BUTP-99/16, Bicocca-FT-99-25 
Effective Lagrangians in bound state calculations 



V. Antonelli a ' fe , A. Gall"*, J. Gasser a , and A. Rusetsky a ' c 
a Institute for Theoretical Physics, University of Bern, Sidlerstrasse 5, CH-3012, Bern, 

Switzerland 

b Dipartimento di Fisica, Universitd di Milano e Milano-Bicocca e INFN, Sezione di Milano, 

Via Celoria 16 - 22100 Milano, Italy 
c HEPI, Tbilisi State University, University St. 9, 380086 Tbilisi, Georgia 

(13 March 2000) 



Abstract 

In order to investigate specific aspects of bound state calculations in a non- 
relativistic framework, we consider the energy-levels of a massive scalar par- 
ticle, which moves in an external field and interacts in addition with a mass- 
less scalar particle. The discussion includes the following topics: dimension- 
ally regularized bound-state calculations, ultraviolet finiteness of bound-state 
observables and their independence of the off-mass-shell behavior of Green 
functions, non-renormalizable interactions, structure of the non-relativistic 
two-point function, power counting and matching. 
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I. INTRODUCTION 



In several recent publications (fTJ-0], the decay of 7r + 7r~ atoms has been studied in the 
framework of a non-relativistic effective Lagrangian - a method originally proposed by 
Caswell and Lepage |4[] to investigate bound states in general. This method has proven 
to be far more efficient for the treatment of loosely bound systems - such as the 7r + 7r~ atom 
- than conventional approaches based on relativistic bound-state equations ||. In JJ, we 
have discussed the chiral expansion of the decay width in this framework. The corresponding 
numerical analysis is carried out in 0. 

While performing that investigation, we have been confronted with the following specific 
features of this approach: 

i) The non-relativistic effective Lagrangian defines a non-renormalizable Quantum Field 
Theory, where Green functions are expanded in powers of the coupling constant and of 
external momenta. While the expansion in the coupling constant can be done in the 
standard manner, the expansion in momenta is highly non-trivial: In order to make 
sense of the perturbative expansion, one needs a procedure where the contributions of 
higher dimensional operators in the loop expansion are suppressed. This problem goes 
under the name "power counting" . 

ii) The renormalization program in this non-relativistic theory must be carried out in a 
consistent manner. A reduction formula should be available, that allows one to evaluate 
scattering matrix elements from the residues of off-shell Green functions, analogous 
to the procedure used in the relativistic formulation. For this, the structure of the 
two-point function in the non-relativistic theory must be known. 

iii) One needs a systematic framework to calculate energy-levels and decay widths in 
the non-relativistic theory. In view of applications in hadronic atom calculations, 
the framework should allow one to also evaluate the lifetime of unstable states in an 
unambiguous manner. 

iv) The effective Lagrangian contains terms with arbitrary high powers of space deriva- 
tives. As a result of this, the matrix elements of the perturbation Hamiltonian between 
the unperturbed wave functions start to diverge at a sufficiently high order of the per- 
turbative expansion. One should demonstrate that these divergences cancel at any 
given order in the expansion. 

v) The scattering sector of tttt interactions is most effectively described by Chiral Per- 
turbation Theory (ChPT). Here, the ultraviolet (UV) divergences are removed by 
higher order counterterms in the chiral expansion. One needs an analogous statement 
concerning the ultraviolet divergences in bound-state observables. 

Some of these issues - like the problem of power counting in the context of QED bound 
states [|HU, NRQCD and HQET g-pf, effective field theories in the one and two- 

nucleon sectors |16| - have been treated previously in the literature, while some others have 
not yet been, as far as we are aware, considered in detail. In particular, a detailed discussion 
of the renormalization program in the non-relativistic sector - including an investigation of 
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the structure of the two-point function - is not yet available to the best of our knowledge. 
For this reason, we present in this article two examples of Quantum Field Theories, where 
these issues can be studied in an unambiguous manner. 

The article is organized as follows. In section |J, we present the relativistic formulation 
of the models: First, we consider a massive scalar particle moving in a given external field. 
Second, we add a Yukawa interaction of the heavy particle with a massless scalar particle. 
This interaction changes the energy-levels - we evaluate the shift at lowest non-trivial order 
in the ground state. This section also contains an example of a non-renormalizable interac- 
tion, that nevertheless results in ultraviolet finite shifts of the energy-levels. In section [TTI|, 
we formulate the non-relativistic version of the external field problem. A general formula for 
evaluating bound-state poles is derived with the use of Feshbach's [|l7j projection technique. 
We determine the perturbative expansion of the energy-levels, using dimensional regular- 
ization to tame the ultraviolet divergences that occur in these calculations. In section [IV], 
we discuss the inclusion of a dynamical light field in the non-relativistic framework, and 
illustrate the problem generated by loop corrections. The section [V] is devoted to a detailed 
study of the two-point function in the non-relativistic theory, while the matching procedure, 
needed to determine the coefficients in the non-relativistic Lagrangian, is considered in sec- 
tion [VI|. In section |V1J] , we evaluate the energy-level shift in the ground state and explicitly 
verify that the result agrees with the relativistic calculation performed in section [H]. Finally, 
section |VII]] contains a summary and concluding remarks. The Appendices are devoted to 
technical details of the calculations presented in the text. 



II. RELATIVISTIC FORMULATION 



We investigate in this section the bound state spectrum of a heavy particle, which is 
bound in an external Coulomb field, and which interacts in addition with a massless scalar 
particle. The results of the relativistic formulation outlined here are compared in later 
sections with a non-relativistic framework. 



A. The external field 



We consider a heavy, complex scalar field H of mass M, interacting with an external, 
static source 0, 

C = d fl Wd^H-M 2 H^H + \H^H(j). (2.1) 

To determine the energy spectrum, one may proceed in several ways. In the following, we 
consider the connected two-point function 

G{x,y)=i(0\TH{x)Hi{y)\0) e , (2.2) 

because its Fourier transform contains the bound state poles. It is useful to note that G is 
the inverse of a differential operator, 

(D M - \<j>)G(x, y) = 5\x - y) ; D M = U x + M 2 , (2.3) 
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result of which one has[] 

(■ I hi ■ Ihi'l'Ih, , 

T = (1 — X^Dlj-y 1 X(j) = \<j) + \ 2 <j)D^(f) + 0(A 3 ) . (2.4) 

The scattering amplitude is obtained from the Fourier transform of T by putting the mo- 
menta on the mass-shell, 

J dxWxrfyj^-^e-^j^Tfay) = T(p, k) , (2.5) 

with 

p 2 = M 2 , |p| = |k| . (2.6) 
Expanding in powers of A, we find 

T(p, k) = £ A n / # n (p, k; q l5 . . . , q n ) J] 0(q,)^(q t ) , (2.7) 

n=l J i=l 

with 

0(q) = y rf 3 xe" iqx 0(x) , du(q) = c/ 3 q/(2vr) 3 . (2.8) 

The kernels K n are independent of the mass M, 

^ 1 (p,k;q 1 ) = (27r) 3 5 3 (q 1 -p + k), 

K 2 (p, k; qi , q 2 ) = (27r) 3 5 3 ( qi + q 2 - p + k) 1 . (2.9) 

(k + q 2 ) 2 - k^ 

In the following, we often consider the case where the external field is of the Coulomb-type, 

0=j^p ( 2 - 10 ) 

Below, we shall not always distinguish between the case of short and long range (Coulomb) 
external sources. In the latter case, the on-shell scattering amplitude contains infrared 
divergences that are tamed by using dimensional regularization. 
The Fourier transform of the propagator 

G(E;p,k) = y c /x°rf 3 xrf 3 ye ii?(:c0 - J/0) e- ip V ky G'(x,i/) (2.11) 



is then closely related to Schwinger's Green function |I8 



1 We use throughout the Landau symbols O(x) [o(x)] for quantities that vanish like x [faster than 
a;] when x tends to zero. Furthermore, it is understood that this holds modulo logarithmic terms, 
i.e. we write also 0{x) for xlxix. 
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G (g;p , kH -i^ Gs ( g2 - M2 ;p , k ) . (2.12) 

2M \ 2M Jz,-*=A, 



' 2.1/ 



The explicit expression is given byQ 

G(A, p, k) — — — — + 



£2 _ M 2 _ p 2 £2 _ M 2 _ p 2 | p _ k |2 £2 _ M 2 _ k 2 

+ V - M 2 - p> toA " /( ^ P - k) g»-^-k» • (2 ' 13) 



with 



I(E;p,k)=J Q [( p _ k )2 p + z/ 2/ A 2( 1 _ p )2( M 2 + p 2_ E 2^ M 2 + \g2 _ E 2^ ' ( 2 " 14 ) 

where z/ = | A (M 2 — E 2 )~ 1 / 2 . 

In the representation ( |2.13| ), the two-point function is split into three parts. The two 
first terms correspond to free propagation and to one interaction with the external field, 
respectively, whereas the third term describes multiple interactions. The function G(E; p, k) 
develops poles at v = n or, equivalently, at 

£ 2 = M 2 -^. (2.15) 

The residues factorize in a product of two wave functions which are solutions of the 
Schrodinger equation with a Coulomb potential, 



E 2 n — M 2 — p 2 Vfe(p) = -4ttA / ^M^ n (q) . (2.16) 

p — q| 



For later use, we display the normalized ground-state wave function in momentum space, 
MP) = \ 2 1 > ,7 = A/2 ; J diz(p) |^(p)| 2 = 1 . (2.17) 



B. Radiative corrections 

We now supplement the Lagrangian fl2.1|) with a Yukawa-interaction of the heavy particle 
with a massless scalar field. This additional interaction affects the energy-levels ( |2.15| ) - we 
shall investigate the induced change in E\ at lowest non-trivial order in the Yukawa coupling. 



2 In order to simplify the notation, we omit the negative imaginary part in M 2 that is understood 
here and in the following. 
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1. The Lagrangian 



The Lagrangian is 

C = d^H^d^H - M 2 H^H + -(d^W/ - m 2 i 2 ) + WH{\<f) + e£) + C c . t . , 

£ ct = -5M 2 H ] H - ^Sm 2 £ 2 + (c + c x A0 + c 2 A 2 2 K • (2.18) 

The mass m is used as an infrared regulator in intermediate steps. Whenever possible, we 
always consider the limiting case m — > 0. The counterterms - collected in the quantity £ c t . 
- cancel the ultraviolet divergences in the diagrams displayed in Fig. |l|. As is discussed 
below, the counterterm proportional to c 2 is needed in order to have a smooth behavior of 
the scattering amplitude in the limit m — > 0. 
In dimensional regularization, we find 

5M 2 = -e 2 Z M + ^- 2 +0(^ 7 ) , 5m 2 = -e 2 L M + 0(^L\ 



eM 2 - e 

Co = -eM 2 L M + —, c 1 = eL M , c 2 = , (2.19) 

where D stands for the dimension of the space-time, and 

1 ( = , , 1, M 2 - 



7T Z \ 2 fi Z 



m = V D ~ 4 - ^r'(l) + ln47r)) , (2.20) 

where /i denotes the renormalization scale. The finite parts of the mass counterterms are 
chosen such that M coincides with the physical mass of the heavy particle at order e 2 , and 
that the light field stays massless at this order. As is discussed below, the finite parts of the 
couplings Q as chosen in (|2.19|) guarantee a smooth behavior of the scattering amplitude 
in the limit m — > 0. With ( |2.19| ), the Green functions of the heavy and light fields are 
ultraviolet finite at order e 2 and at any order in A. 



2. Perturbation theory for bound states 



The Yukawa interaction ( |2.18| ) generates a shift in the energy-levels ( [2.15| ). In order 



to derive the formula that allows us to calculate this shift, we start from the connected 



two-point function of the heavy scalar field in ( 2.18|) , 

G(x,y) = z(0\TH(x)H\y)\0) c . (2.21) 
It obeys the familiar Dyson- Schwinger equation 

G{x,y) = G(x,y)+ f dudvG{x,u)E(u,v)G(v,y) , (2.22) 
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where G stands for the Green function ( |2.2|) in the absence of the Yukawa interaction. The 
relation ( j2.22j ) defines the self-energy E of the heavy scalar particle. The expansion of E in 
powers of the coupling constant e takes the form 

E(w, v) = £ 2 (u, v) + E 4 (w, v) + --- , (2.23) 

where E n = 0(e n ). Explicitly, we find for the lowest order expression 

E 2 ("U, v) — e5(u — v) J dwD^{u — w) [— ieG(w, w) + c + C\\<j)(w) + c 2 \ 2 (fi 2 (w)] 

- 5M 2 5(u-v) -ie 2 D m 1 (u-v)G(u,v) , (2.24) 

where D^-(x) denotes the propagator of the free two-point function of the light field, 
(a + m 2 )D m 1 (x) = 5\x). 

It is convenient to define the Fourier transforms 

E(£; x, y) = J dx Q e lE{x ° ~ y0) E(x, y) , (2.25) 

and 

E(E; p, k) = J dVy e-* px+iky E(£; x, y) . (2.26) 

The diagrams contributing to E 2 are depicted in Fig. 0. We have not displayed the 
contributions generated by the counterterms Cj and SM 2 . There are two types of diagrams: 
the tadpole-like diagram I, where the loop is generated by the propagator G, and the self- 
energy diagram II, generated by a closed loop with the light scalar and G. The tadpole-like 
diagram leads to a self-energy of the form 

E^E; p, k) = e 2 -^ ; q = p - k , (2.27) 

where the residue behaves at low momentum transfer as 

R = R 5 3 (q) + ^ + ^ + R 3 ln^ + O(l) , q^O , (2.28) 
q |q| M 

and where R m = 0(X m ). The finite parts in the counterterms are chosen such that -Ro,i,2 = 0. 
As a result of this, the regulator mass m can be sent to zero, and the self-energy is not more 
singular (except from logarithms) than the external Coulomb-field in the forward region 
p = k. In particular, one has 

P ' k) = "24^4 f 1 + ° (q2) } + ° (e2A2) • (2 ' 29) 
The coupling to the external Coulomb-field is modified accordingly at this order, 

„2 



A -> A 



1 - 



96tt 2 M 2 



(2.30) 
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Note that in QED, the coefficients Ro,i,2 are identically zero by charge conjugation and gauge 
invariance. 

Coming back to the bound-state energy, we observe that the Fourier transform of G(x, y) 
has a pole at the energy E B corresponding to the bound state B, 



G(£;x,y) = J dx°e iE ^~^ G(x, y) 



(0\H(0,x)\B)(B\W(0,y)\0) 
2E B {E - E B ) 



+ ■ 



(2.31) 



Using Eqs. (|2.3|), (|2.22| ) and ( |2.31| ), one obtains for the bound-state wave function, defined 
as VK X ) = (0|JT(0, x)|-B), the equation 



E 2 B -M 2 + A X + A0(x) ) ^(x) = ~Jd 3 y Z(E B ; x, y) ^(y) . 



^2.32) 



In order to calculate the corrections to the energy-levels, one can apply ordinary pertur- 
bation theory to the Schrodinger-like equation ( f2.32j ). The result is 



K-K 



dv{p) dv(k) C(P) ^2{E n] p, k) i> n (k) + 0(e 4 ) . 



(2.33) 



Here ipn denotes the normalized unperturbed Coulomb eigenfunction ( j2.16| ) , whereas E 2 and 
E 2 stand for the true and unperturbed eigenvalues, respectively fL9f . 



3. The ground state 

We are now ready to evaluate the shift in the energy- levels, generated by the Yukawa 
interaction eWHl in (|2.18|) . We restrict ourselves to the calculation of the shift in the 



ground-state. By dimensional reasons, one has 

AE 2 = e 2 F(z) + 0(e 4 ) , z = X/M . (2.34) 

Below, we determine the leading term in the expansion of F in powers of z. The calculation 
can be performed at zero mass of the light particle, because the bound particle is always 
slightly off-shell - there is no additional infrared regulator needed. The shift is naturally 
split into the two parts originating from the diagrams of type I and II displayed in Fig. |2|. 
We split the contributions from the diagrams of type II furthermore into a zero-Coulomb, a 
one-Coulomb and a multi-Coulomb part according to Eq. ( p.!3| ), and write 

F = Fj + Fj/,0 + F u>1 + F n , m • (2.35) 

The contribution from the diagrams of the type I is remarkably simple to evaluate. Indeed, 
the previous discussion has shown that - at leading order in A - the inclusion of these 
diagrams merely amounts to an 0(e 2 ) correction in the strength of the Coulomb potential, 
see Eq. ( |2.30|) . The contributions of order A 2 and higher in £2 ma y be discarded at leading 
order. One therefore has 

z 2 

F7 = + o(z 2 ) . (2.36) 
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The expression for the zero-Coulomb part in the diagram II starts again at order z 2 , mod- 
ulated with a logarithmic term, 



z 2 / I , 



The one-Coulomb contribution is given by 

Fu,i = - / du(p)du(k) VI(p) r(£7 x ; p, k) ^(k) , (2.38) 

where T stands for the one-loop correction to the coupling of the heavy particle to the 
external field, 



r(£ i;P ,k) = - J 



d*q 1 1 



(2tt)H q 2 (M 2 - {Ex - q ) 2 + (p - q) 2 )(M 2 - (E 1 - q ) 2 + (k - q) 2 ) ' 

(2.39) 



The calculation of the one-Coulomb contribution to the energy shift is relegated to Ap- 
pendix [A]. We simply quote the final result, 

z 2 

F n ,i = — ^ (In z 2 - 6 + 8 In 2) + o{z 2 ) . (2.40) 

The logarithmic terms cancel in the sum of the zero- and one-Coulomb contributions. 

The calculation of the multi-Coulomb contribution is given in Appendix |TJ. The result 
can be written as 



z 2 



II, m 



327T 5 

with 



J + o{z 2 ) , (2.41) 



J = 2567T 2 / r-r-7- / du(v)dv(v) / r-TTT ^-77 — tt X 

Jo (l + xy/ 2 J K J V 'Jo (l + u 2 ) 2 (l + v 2 ) 2 (l + x + u 2 )(l + x + v 2 ) 

1 (2.42) 



(u - v) 2 p + (1 - p) 2 (l + x + u 2 )(l + x + v 2 )/(4(l + x)) 
Collecting everything, we find 



F(z) = -^(y-8 ln2 + J) + o(z 2 ) . (2.43) 



C. Non-renormalizable interactions 

According to Eq. ( |2.33| ), the corrections to the Coulomb energy-spectrum are obtained 
at leading order by folding the self-energy with the unperturbed Coulomb wave functions 
in momentum space. In theories with non-renormalizable interactions - e.g. in ChPT - the 
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amplitudes grow faster in the external momenta at every order in the loop expansion, and 
the just mentioned integration starts to diverge at some order. This divergence is not, at first 
glance, cured by the counterterms in the Lagrangian, because it arises from an additional 
integration of the matrix element that was already made finite by these counterterms. In the 
context of the 7r + 7r~-atom, it has been claimed in Ref. |2D) that, for this reason, one needs to 
introduce an explicit ultraviolet cutoff in order to render the bound-state observables finite. 
Here we wish to demonstrate in a simple model - bearing all the important properties of 
the original problem - that this is not the case: the ultraviolet divergences that occur at 
intermediate steps in the evaluation of the perturbed energy-levels are all cancelled. 
In order to illustrate the point, we consider the Lagrangian 

C = d^d^H - M 2 H ] H + \H ] H(j) + g {(□ + M 2 )H^} A0 {(□ + M 2 )H} + C c . t . , (2.44) 

where is a Coulomb field. In the last term, we anticipate the necessity of adding countert- 
erms to cancel the ultraviolet divergences in the Green functions. 

The self-energy part is again defined by ( |2.22|) , and the shift in the energy- levels can be 



calculated from (|2.33| ), with e 2 —>■ g. At order g, we find for the self-energy part 

E(£; p, k) = (M 2 - E 2 + p 2 ) (M 2 - E 2 + k 2 ) + £ c ,.(£; p, k) . (2.45) 

IP k l 

The first term is ultraviolet finite - however, it leads to an ultraviolet infinite contribution 
to the energy shift at order g. This fact mimics the situation which one encounters in the 
calculation of the energy of hadronic atoms beyond the so-called "local" approximation |20 . 
The remedy for the difficulty is provided by the counterterm in Eq. fl2.44| ): as one may 
observe from Eq. ( |2.22|) , the two-point function [defined with the Lagrangian ( |2.44|) 1 is infinite 



at order g, although the self-energy is not. In dimensional regularization, the ultraviolet 
divergent piece in the two-point function is given by 



G(E;p,k) 



div 



-2n\ 3 gL(n) J dv(q)dv(\)G(E; p, q) G(E; 1, k) + 0(g 2 ) , (2.46) 
where L(fj,) is given by 

L( yU ) = ( / i 2 ) B - 4 f-i--r(l)-ln47r") . (2.47) 



,£>-4 

This divergence is cancelled by the counterterm 

£ c .t. = ~\ 3 gcH^x)H(x)A<j ) (xy, c = L(jjl) + c r (^) , (2.48) 

where c r (/x) denotes the arbitrary scale-dependent finite part of c. If one would carry out 
the perturbative expansion of the two-point function to higher orders in g, one would have 
to introduce additional counterterms, containing higher derivatives - the situation clearly 
resembles the one in ChPT. Evaluating S ci . generated by the counterterm ( |2.48 ), one can 



verify that the corresponding contribution to the energy-level shift cancels the ultraviolet 
divergence generated by the first term in (|2.45| ), leading to a finite and scale-independent 
result at order g. 
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Let us note that the ultraviolet divergences can be removed from the bound-state observ- 
ables in a perturbative manner: indeed, the divergent piece of the two-point function (p. 46 ) 
was calculated at 0(g). The corresponding counterterm removes the ultraviolet divergence 
in the energy- levels at 0(g). This suggest that also in non-renormalizable theories, bound- 
state observables are rendered finite at a given order, provided that the renormalization in 
the scattering sector is carried out at the same order. 



III. NON-RELATIVISTIC FORMULATION: 
THE EXTERNAL FIELD 

In the remaining part of this article, we wish to reformulate the relativistic problem 
considered in section |I| in a non-relativistic framework. The energy-levels can then be 
calculated using Rayleigh-Schrodinger perturbation theory. At the same time, the formalism 
developed allows one to calculate the width of decaying states in a crystal clear manner, 
without using handwaving arguments. The procedure goes in two steps ||]. First, one 
constructs an effective theory that allows one to evaluate the scattering matrix elements in 
the low-energy region. Second, one uses the Hamiltonian of this effective theory to evaluate 
the bound state poles. 

We restrict the considerations in this section to the non-relativistic formulation of the 
external field problem (|2.1|) . Radiative corrections generated by the Yukawa interaction 
( [2.18| ) will be treated afterwards. 



A. The scattering sector 

For the construction of the non-relativistic effective Lagrangian that is relevant for the 
theory defined in Eq. (|2.1| ), one may proceed in several ways. Below, we choose a method 
that is based on the following observation. In the relativistic framework, the scattering 
matrix elements are obtained from the residue of the poles in the Fourier transform of the 
two-point function G. These poles are generated, in the present case, by the free propagators, 

d m(p) = Ti^—i = ^- (^-^ + — ^1 ; ^ = ( m2 + p 2 ) 1/2 > (s- 1 ) 

which have poles at p° = ±u; p . It is therefore sufficient to construct a non-relativistic theory 
that has the same structure in the vicinity of p° = u p . The pole terms in ( |3.1|) can be 
generated with the inverse of the differential operators 



D ± = ±id t - Jm 2 - A , D ± (x)A ± (x) = 5 {4) ( 



x 



f d^p c ^ x 

A±(X) = J J2^ ±p«-VM 2 + P 2 ■ (3 ' 2) 

The Lagrangian will furthermore be quadratic in the non-relativistic heavy field h. Writing 
£-nr = h^Vh, the differential operator must V contain, in addition to D + , further terms 
that generate the transition amplitude T in ( |2.4|) . Writing 
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V = D + + X 



V- 1 = D- 1 - D- + l T NR D- + l 



(3.3) 



and requiring that T NR is closely related to the relativistic scattering amplitude T in 
one finds [21] that one may take 



X = Xd(j)d(l + \DZ 1 d(fid)~ 1 



d = (2a/M 2 - A)~ 1/2 



with 



Ti 



NR 



dTd . 



(3.4) 



(3.5) 



This formulation involves a non-local Lagrangian that generates the same truncated off- 
shell two-point function as the relativistic theory. In order to arrive at a local Lagrangian, 
one may expand the non-local operators Dg 1 and d in inverse powers of the heavy mass 
M. In this manner, one arrives at a local Lagrangian, that contains, however, an infinite 
string of terms. Furthermore, it contains time derivative of arbitrary high orders, acting on 
the heavy field. These time derivatives may be eliminated in a standard manner by use of 
the equation of motion (EOM). We find that this can be shown in a transparent manner 
by use of external sources coupled to the heavy field. We relegate these manipulations to 
Appendix y and simply quote the result, 



C NR = tfD Q h+Y i 



n=0 (2M) 
D = id t -M+ A/2M 



2n+l 



c 



2n+l j 



C x = rf\<j)h , C 3 = h) A 2 + X((f)A + A0) + \ 2 <f) 2 



C 5 = -h) [4A 3 + A (2A0A + 50A 2 + 5A 2 



3 j.3 



h 



h . 



(3.6) 



+ h ] [A 2 (20 2 A+A0 2 + 30A0) + 2A 

The scattering matrix elements are obtained by evaluating the two-point function 

G NR (x,y) = i(0\Th(x)tf(y)\0) (3.7) 

in the presence of the interactions £j. The perturbation theory is performed as follows. 
First one sums up all mass insertions in the external lines. The Fourier transform of the 
two-point function then develops poles at p° = u> p and at p° = uj-^, 



G NR (p°;p,k) = J dx°d 3 xfi 3 ye ip0(a:0 ^ 0) e- ipx e iky G' J v jR (x,y) 
(2tt) 3 5 3 ( P - k) 



1 



1 



T NR (p,k;p 



1 



1 



(3.8) 



Up -p° LU p - p° y^ZuJ^ V2U\ LU k - p° 

With the normalization chosen, T^ R coincides with the low-energy expansion of the rela- 
tivistic amplitude T in (2.5) order by order. The presence of the factors y/'2ujp, y/2oJ^ in the 
above formula reflects the relation (|3.5|) . The scattering amplitude is obtained by putting 
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the momenta on the mass-shell p° = u p , |p| = |k|. Expanding in powers of the coupling 
constant, one has 



oo „ n 

T NR (p, k) = ^ A™ / fc n (p, k; qx, . . . , q„; M) J[ 0(q;)c?y(q 

n=l J i=l 



(3.9) 



In the internal lines, one need not fully sum the mass insertions. Instead, one sums, at a 
fixed order in A, all terms at a given order in 1/M, 



1 



kn(p, k; qi, . . . , q„; M) = J2 j^ k n,v(P> k; qi, • • • , q 



u=0 



(3.10) 



In Fig. |] we display the corresponding graphs at order A 2 . The diagram Fig. |HA stands 
for the Green function evaluated in the relativistic theory, whereas mass insertions and the 
contact term in the non-relativistic theory are displayed in graphs Fig. §B and Fig. |3]C. The 
mass dependence cancels, order by order, and the leading term coincides with the matrix 
element in (|2.9| ), evaluated in the relativistic theory, 



k 



n,0 



K ■ h ■ 
- I1 n ; i^riji 



,i > 1 . 



(3.11) 



We illustrate the procedure with the matrix element at order A. The first two terms in the 
1/M expansion give 



T^(p,k)=A0(p-k)^ 



2M 2 



0(M~ 



+ 0(A 2 ) . 



(3.12) 



A different Lagrangian would be generated by use of the standard Foldy-Wouthuysen-Tani 
transformation, which gives 



FWT 



h\id t 



M 2 - A - \<f)] 1/2 )h . 



(3.13) 



Again, the Lagrangian may be expanded in inverse powers of M. At order 1/M 5 and higher, 
it differs from the one considered above, 



C 



NR 



A 



FWT 



QAM 5 



h) { A 2 (f) - 2A0A + 0A 2 - 2AA0 2 + 2A0A0} h + 



(3.14) 



On the other hand, evaluating the scattering amplitude with (|3.13| ), one finds that it agrees 
with the previous one on the mass-shell, to all orders in the coupling constant. The reason is 
the fact that the two Lagrangians may be transformed into each other by a field redefinition. 
To illustrate, the terms at order 1/M 5 are removed if in Cnr we substitute 



h (1+X)h 

x- 1 



32M 4 



A(A0 - <j)A) - A 2 2 



(3.15) 



On the other hand, the two-point functions differ off the mass-shell. Indeed, consider the 
term linear in A. From ( |3.14| ) it follows that the difference is proportional to 
X(j)(p - k)(p 2 - k 2 ) 2 at order M" 5 . 
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B. Perturbation theory for bound states 



In Quantum Mechanics, bound states are often treated by splitting the Hamilton operator 
into a part with known spectrum and known wave functions, and a perturbation whose effect 
can be systematically calculated in a power series in some small parameter. This procedure 
is successful, because all truly non-perturbative effects are collected in a compact manner 
in the lowest order wave functions. In Quantum Field Theory, however, the basic quantities 
are the Green functions, and there is no simple analogue to wave functions. This makes 
the treatment of bound states inherently more complicated. The salient feature of the 
effective non-relativistic theory is that it allows for a consistent Hamiltonian formulation in 
the context of a field theory. 

To see how this works, we define the Hamiltonian pertaining to the Lagrangian (|3.6| ), 

Ti = Tio + Tic + 7~ti, 

The static external field preserves the time invariance of the theory. We perform canonical 
quantization at t — 0, introducing creation and annihilation operators 

a t(p) = J d 3 xe ip *rf(0,x), a(p) = / rf 3 xe" ipx /i(0, x) , (3.17) 

obeying the commutation relation [a(p), a^(q)] = (27r) 3 5 ( - 3 ^(p — q). The Hamilton operator 
is defined by H = J c/ 3 x7i(0,x) and similarly for Ho, He and Hi. A free non-relativistic 
particle is described by the state |p) = a^(p)|0), which obeys H |p) = (M+p 2 /2M)|p). The 
Coulomb wave functions ip n are defined as the solutions of the non-relativistic Schrodinger 
equation, 

p 2 \ . , , 2ttA 



M+ 2M) MP) ~~WI du ^]^T~^ Mq) = ^»<P)' (3- 18 ) 



with 




En = M-— — (3.19) 



for the discrete spectrum. Note that the wave functions ( |2.16| ) satisfy the equation ( |3.18| ). 
The eigenstates (Ho + Hc)|n) = E n \n) of the unperturbed Hamilton operator are given 

by 

\n)= Jdu(p)Mp)\p)- (3-20) 



Our task is to find the isolated poles of the two-point function ( |3.8| ). Here we note that 
the spectrum of the full theory is also encoded in the singularities of the resolvent 

G(E) = ^-Lj . (3.21) 
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If the Hamilton operator annihilates the vacuum, H|0) = 0, the resolvent is closely related 
to the two-point function, 

(p\G(E)\k} = -G NR (E;p,k) , (3.22) 

see Appendix |D|. The advantage of using the resolvent instead of the Green function is due 
to the fact that the bound state poles can be isolated in a very convenient manner. As a 
lowest order approximation, we ignore Hj. The resolvent 

G(E) = 1 (3.23) 

hi — II — ±lc 

of the unperturbed Hamilton operator is proportional to the Coulomb Green function given 
by Schwinger's solution, 



(p|G(£)|k) = (27r) 3 G 5 (£-M;p,k) 



^ e 2M 



(3.24) 



If the perturbation is small, the poles of G(E) are expected to lie in the vicinity of those 
of G(E). To work out the exact position of the poles, it is convenient to first consider the 
operator Ti defined through 

Ti(£) = H : + HiG(£)Ti(£) . (3.25) 

It is related to the full resolvent by G(E) = G(E) + G(E)T l (E)G(E). Applying a technique 
introduced a long time ago by Feshbach |nj, we use the projector \n)(n\ to exhibit the pole 
near the unperturbed value E n , 

T ,„ .„ , m , T„(E)\n)(n\T n (E) 

Tl(B) = T " {E) + E-E n -(n\T n (E)\ny ^ 

where T n is defined by 

T n (E) = Hj + HiG n (E)T n (E) , G n (E) = G(E)(1 - \n)(n\). (3.27) 

The corresponding decomposition of the full resolvent is obtained by expressing G(E) 
through T n (E), 

G(E) = G n (E) + G n (E)T n (E)G n (E) + (1 + g^lglM^ -± ^ T 
v ; y ' n K j uk j n K j E - E n - (n\T n (E)\n) 

(3.28) 

The energy of the bound state is determined by the pole position E n in the full resolvent. 
Since both G n (E) and T n (E) are regular in the vicinity of E n that includes also E n , the 
pole position is given by the zero of the denominator in the last term of Eq. (|3.28j) [ |17|j , 

E n — E n — (n\T n (E n )\n) = 0. (3.29) 

So far, we did not make any approximations or pose restrictions on the perturbation Hj. In 
fact the last formula holds even when we add a "radiation field" , see Section |VII] . Moreover, 
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the same formula can be used for the analysis of the decays of quasistable bound states 
- without introducing the concept of "wave function of the unstable state" . The problem of 
determining the positions of the poles in the two-point function Gnr(E; p, k) is thus reduced 
to the calculation of the matrix elements (n\T n (E)\n) in old-fashioned perturbation theory. 

Having taken the non-perturbative part into account, it remains to solve Eq. ( |3.29| ) in 
perturbation theory. In the case at hand, A/2M is the only small parameter and we can 
express E n as a power series 

oo 

E n = E n + ^ ^ m Cmn- (3.30) 

m=3 

The coefficients C mn can be obtained by using in ( |3.29|) the iteration 

(n\T n (E n )\n) = (n^n) + (n|H I G n (^ n )H I |n) + • • • , (3.31) 

which can be calculated by inserting a complete set of states - corresponding to H + H c - 
and keeping terms of the appropriate order in A. Here, we recognize the familiar Rayleigh- 
Schrodinger perturbation theory. 



C. Energy- levels 

With a given Lagrangian, and a framework set, one can now calculate the energy of the 
bound-state levels in the non-relativistic theory - the final answer should reproduce of course 
the relativistic result. For example, for the ground-state energy in the non-relativistic theory 
one should obtain order by order in perturbation theory in the coupling constant A 

El= M 2 --) =M-- ^__^_ + .... (3.32) 

V 4 7 8M 128M 3 1024M 5 v ; 



Using Eq. ( [3.31|) , one encounters two conceptual difficulties. First, as was demonstrated 



in the preceding discussion, the non-relativistic Lagrangians which can be derived from the 
relativistic theory, are not unique, although all of them yield the same S'-matrix elements. 
In the case considered here, they differ by terms which vanish by the use of EOM and field 
redefinitions. Second, although all couplings in the non-relativistic Lagrangian are finite, it 
contains terms with arbitrary high powers of space derivatives. This means that the matrix 
elements of the perturbation Hamilton operator between the unperturbed wave functions 
start to diverge at a sufficiently high order of the perturbative expansion. Consequently, 
one has to demonstrate that: i) the terms which vanish by the use of EOM or which 
can be removed by field redefinitions, do not contribute to the bound-state energy, and 
ii) divergences in the bound-state energy cancel when all terms at a given order in the 
coupling constant are summed up. We shall ensure that this is what indeed occurs in 
explicit calculations. 

In order to evaluate the energy-level shift, one has to specify the explicit form of the 
Hamiltonian. We find it useful to start with the FWT Hamiltonian, which has a remarkably 
compact form, 

n FWT = V^-A-A^ h = h^M- ^t)h - (^g^)! + ...),. (3.33) 
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The first two terms in Hfwt = Hq + Hc + 'Hi correspond to the unperturbed non-relativistic 
Coulomb problem, and the rest is considered as a perturbation. In order to regularize 
the ultraviolet divergences mentioned, we perform the calculations in d space dimensions 
[d = D — 1], and perform the limit d — > 3 at the end. We define the Hamilton operator 
Hfwt = / d d x.H.FWT(0, x), and similarly for H , He and Hi, etc. Further, in d dimensions, 
the momentum-space Coulomb wave function ip d (p) is chosen to obey the <i-dimensional 
Schrodinger equation with the same kernel as in Eq. ( |3.18| ), 

[K - M - t.) Hip) = ~ j dMl) i^ji rf(q> , (3.34) 

where di/^q) = d d q/(2n) d , and E d denotes the energy of the n-th eigenstate in d dimen- 
sions^ Using Eq. (|3.34|) , it is seen that the unperturbed Coulomb wave functions are eigen- 
functions of Hi, 

J d d x/i t (A + \(f)) k h\n) d = (2M(M - E d )) k \n) d , (3.35) 

where \n) d is constructed with the use of Coulomb wave function ^j(p), similarly to 
Eq. ( ggg ). In this expression, one can safely put d = 3. Using then the perturbation 



series for the energy-levels, one ensures that Eq. ( |3.32| ) holds in all orders in the coupling 
constant A. 

The simplicity of the above derivation is a consequence of the choice of the Hamiltonian 
in the FWT form ( |3.33|) and should not be misleading. As we have already mentioned, in 



general one encounters ultraviolet divergences in individual matrix elements in sufficiently 
high orders of the non-relativistic expansion. To illustrate this fact, we again consider the 
original Lagrangian ( |3.(jp . The difference 57ii between the pertinent Hamiltonian ( |3.(j| ) and 
the FWT Hamiltonian ( |3.33| ) contains terms at 0(M~ 5 ) and higher, 

SHj = h\<pA 2 + A 2 - 2A0A)/* + ^ />t(A0 2 - 0A0)/> + 0(M~ 7 ) . (3.36) 

It is expected that the level shift induced by the perturbation blti vanishes to all orders 
in A. Consider matrix elements of the form d(n\ J d d ^h\/\) k ( y \(j)) l ~ k h\m)d 1 with all possible 
permutations of the operators. Using dimensional arguments, one finds that all such oper- 
ators will contribute at order 0(A 2Z ) to the bound-state energy (up to possible logarithms). 
Therefore, the energy shift induced by STii starts at order A 6 M~ 5 - this term is given by 
the diagonal matrix element of this operator between the unperturbed wave functions. The 
following matrix element is then e.g. needed, 

J d d x d {n\rf A0A/i|m) d = J du d (p) du d (k) ^(p)*p 2 ^ k 2 ^(k) 



-27rA 2 ^(0)^(0)(L(/i) + ln^)+O(l), rf^3, (3.37) 



3 We have checked numerically that the ground state energy depends smoothly on d in the vicinity 
of physical space dimension d = 3. In particular, no bifurcations occur. 
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where 



(3.38) 



and where is given by Eq. (|2.47|) . Eq. (p . 3 7|) illustrates the ultraviolet divergences 

mentioned. Using the bound-state equation in d dimensions ( |3.34| ) in the calculation of all 
matrix elements from ( |3.36| ), one finds that all divergent terms containing L(fi) cancel at 
order A 6 M -5 . Moreover, it can be checked that the diagonal matrix element and, therefore, 
the energy-level shift due to SHi vanishes at this order, 



(n|5Hi|n) d = 0(A 8 M" 7 ) 



(3.39) 



Consequently, at this order the theories described by the Lagrangian ( |3.6| ) and the FWT 
Lagrangian, yield the same bound-state energy. We do not intend to construct a proof that 
this holds to any order. 



IV. INCLUDING THE DYNAMICAL LIGHT FIELD 



In section |HI| , we have set up a non-relativistic framework that allows one to calculate 
the energy-levels of a particle moving in a given external field. We wish to extend this 
construction to the model discussed in section [Tj], and to consider radiative corrections to 
the energy-levels, generated by the Yukawa interaction with the light field. In order to set 
up a non-relativistic framework for the calculation of the energy shift, one has to specify the 
Lagrangian which will be used in these calculations, and the rules that are used to obtain 
observable quantities. Indeed, the perturbative expansion will be carried out not only in 
the coupling constant, but also in energies and momenta, because one needs to know the 
amplitude just in the vicinity of the threshold. In such a multiple expansion, the question 
of ordering arises - one has to specify the relative magnitudes of the expansion parameters. 
This then determines what terms in the expansion should be grouped together. A particular 
ordering is called power counting scheme. The calculational framework must then ensure 
that the ordering is not destroyed by loop corrections - this problem goes under the name 
"validity of power counting" . 

In this section, we discuss scattering amplitudes at tree level and identify the problems 
that show up in loop calculations. In order to keep the arguments simple, we put A = in 
this and in the following section. 



A. Tree graphs 

In a first step, we construct a non-relativistic Lagrangian that reproduces the tree graphs 
generated by the relativistic Lagrangian Eq. (|2.18 ). We will restrict the discussion to the 



one heavy particle sector, and first consider scattering matrix elements in the relativistic 
theory, 

(p, q n+ i, . . . , q m out \k, q u . . . q n m) = (. . .in | . . . in) + z(2tt) 4 5 4 (P / - Pj)T(p, k;q x ,..., q m ) ; 

m>2,n>l. (4.1) 
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At tree level, the matrix element T depicted in Fig. f|, is a sum of products of free propagators 
- each term is of the form 



m—l 



with 

m m 

pf = k° + J2 c v ,d , Pi = k + £ c^q, ; = ±1,0. (4.3) 
i/=l v=i 

The low-energy expansion is now defined as follows. We decompose each propagator in a 
singular and a regular part, 



1 ( 1 1 \ 1 



M 2 - Pf \ Uj - P? ' uij + P° ; 2^- ' Uj ' " ' ' ' 
write in the first part 

p2 

^ = M+ 2M + A ^' = ^ + A ^' ' (4 ' 5) 
and note that Aujj is much smaller than the difference ouj — Pf, because 

Au, = O(PJ) , -P° = 0{$) , (4.6) 

everywhere in momentum space except the specific configurations of external momenta for 
which the difference ujj — P® is of order P* (these configurations have zero measure in 
momentum space). Therefore, one may expand in the first term in Auj. Furthermore, in 
the second term and in the factor 1/uij, one expands in inverse powers of the heavy mass 
M. In this manner, each of the propagators in each factor in the product ( f4.2|) becomes a 
series of terms of non-relativistic propagators. The leading term reads 

{2M) m - 1 ^ a, - Pf ' ^ 

whereas the next terms contain additional powers Aujj/ (ujj — Pf). There are also terms that 
do not contain any propagators - these originate from the second term in Eq. (|4.4|) . Finally, 
we also expand the quantities Aujj in inverse powers of the heavy mass. As an example of 
this procedure, we consider Compton scattering, which amounts to m = 2, n = 1 in ( [4.11 ). 
The intermediate particle in the s-channel graph has four momentum 

If = (* + g?,k+cu) , (4.8) 

and the matrix element becomes 



(il^ + P 2 ) 1 / 2 , (4.4) 



e 2 T, 



2ui 



1 1 

+ 



Ui - k° - q\ ui + k° + ql 



(M 2 + P\) 1 ' 2 . (4.9) 



Expanding in the manner described above, we obtain an infinite series of terms, 
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11 p 2 1 p 4 1 1 
T = — - — + _L! i + _JL_ + . . . (no) 

s 2Miu 1 -k°~q° 1 4M 3 a;i - fc° - 16M 4 (c^i - - g?) 2 4M 2 ' K ' 1 

These contributions have the following representations as Feynman graphs: non-relativistic 
propagation, vertex correction, mass insertion, contact term, etc. The contribution in the 
w-channel may be treated in an analogous way. In this manner, one generates a well-defined 
series of terms, that may be ordered in inverse powers of the heavy mass M. 

The next step consists in the construction of a non-relativistic Quantum Field Theory 
that generates exactly this series for any tree graph in the one heavy particle sector. We start 
from the observation that the tree graphs are generated by the classical action, evaluated at 
the solution to the classical EOM. Retaining only the terms that correspond to the sector 
with one heavy particle, the corresponding non-relativistic Lagrangian is closely related to 
the one discussed in section |TTT| for the motion in the presence of the external field, 

C = i d^t + tfVh , V = D + + ed£d(l + eDZ 1 d£d)- 1 . (4.11) 

The differential operators D± and d are defined in section |m[ In order to arrive at a local 
Lagrangian, one again expands these differential operators in inverse powers of M - the result 
is given by Eq. Q in Appendix 0, with -> el. This La grangian is analogous to the 
effective Lagrangian ( |3.6|) , except that now, the field £ is time dependent. As a result, some 
of the operators contain time derivatives of the light field £. Higher order time derivatives 
can be eliminated in a standard manner, see Appendix |C[ The final result for the effective 
Lagrangian is of the form 

1 OO 1 

C NR = tfDoh + -d,£dH + £ c n , 

d = ehUh, C 2 = 0, £ 3 = rf{A 2 + e{£A + A£) + e 2 £ 2 )h, .... (4.12) 

The scattering matrix elements are obtained by evaluating - at tree level - the connected 
Green functions 

G NR {x, y;z 1 ,...,z m ) = i(0\Th{x)h\y)£{ Zl ) . . . £{z m )\0) c . (4.13) 

The perturbation theory is performed as follows. First, one sums up all mass insertions in 
the external heavy lines. The Fourier transform of Gnr then develops poles, and one ends 
up with a reduction formula for the scattering amplitude that is similar to Eq. (|3.8|) , except 
that there are also poles corresponding to the light fields. This procedure generates the 
low-energy expansion of the relativistic scattering amplitude discussed above. 



B. Loops 

The evaluation of tree diagrams is performed by expanding in inverse powers of the heavy 
mass M. Each term in the Lagrangian contributes in a definite order - the expansion makes 
perfect sense. By calculating loops, one again needs an expansion where the contributions 
from higher dimensional operators are suppressed in a well defined manner - otherwise, the 
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framework is not useful. This problem has obtained an extensive coverage in the literature, 



in the context of different physical problems [pTO , pT2 , 15|| . We first illustrate the problem 



in the case of the correlator of two heavy fields. In the following section, we show that the 
above Lagrangian does lead to a consistent framework, provided that the rules to evaluate 
the Feynman diagrams are properly adapted. 

In order to avoid inessential technical complications, we consider here the two-point 
function at vanishing 3-momentum p in D dimensions, 

G NR (a)=t [d D xe^° (Q\Th(x)tf(Q)\0) = — Q 1 - ; a = 2(M -p°)/M . (4.14) 

J M — p" — E NR (a) 

Here, we have introduced the dimensionless variable a for later convenience, and use the 
quantity T, NR to collect the contributions from the loops, as usual. The dependence of 
the propagator and of the self-energy part on the variables e 2 and M is suppressed to ease 
notation. In the following, we concentrate on the evaluation of Eatr at lowest non-trivial 
order e 2 . Some of the diagrams that contribute at this order are depicted in Fig. |. The 
contribution from diagram Fig. [5]A reads 

^ _^ 1 fd°l 1 



4M 2 i(2n) D J -PM-p° + l + ^' 

Without performing any explicit calculation, we can make the following observations. First, 
the integral is ultraviolet divergent - one has to add counterterms to the original Lagrangian 
to render the propagator finite. The divergence is not suppressed in the large M expansion 
- the counterterm required contributes to the mass term. Further, expanding the integrand 
in powers of a, it is seen that there will be terms analytic at a = 0, as well as a non-analytic 
term ~ a In a, reflecting the fact that the propagator generates a branch point at p° = M. 
The contributions from the diagrams Fig. |5|B + |5|C are given by the integral 

p 2 1 r rl D I 1 
yB+c _ e 1 / 12 f (A 1fi x 

^ NR 8MH(2ir)°J -P M- P o + /o + ^ " ^ ibj 



Again, the divergence of this integral is not suppressed in the momentum expansion - the 
mass term is affected also here. Further, there are contributions analytic in a, as well as a 
non-analytic term ~ a 3 In a. We conclude from these two examples that perturbation theory, 
based on the Lagrangian fl4.12p and using standard dimensional regularization, generates 



a complicated series of terms: in general, all operators in the effective Lagrangian will 
contribute e.g. to the mass. This feature of the standard expansion is very similar to the 
problems already discussed in application to different physical problems [njl2-15 . 



The reason for the breakdown of the counting rules lies in the inconsistency of the 
procedure: the non-relativistic effective theory we have constructed is supposed to be valid 
at very small momenta. On the other hand, using dimensional regularization, we have 
extended the momentum integration to infinity. In the presence of the hard mass scale M, 
the propagators are no longer homogeneous functions of momenta, and the procedure fails 
to suppress contributions from higher dimensional operators. The problem is clearly caused 
by the contributions from large momenta in the integral. Indeed, at very small momenta, 
the term 1 2 /(2M) in the heavy particle propagator may be dropped - then,the integrand is 
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homogeneous in the momenta, and the contributions from higher dimensional operators are 
suppressed, being proportional to <r D ~ 3 (er D_1 ) for the graph depicted in Fig. ||A (Fig. ||B, 
Fig. |5]C). On the other hand, the high energy contributions generate a contribution that is 
analytic at o = 0. One needs a systematic procedure to subtract these unphysical high- mom- 
entum contributions from the loop diagrams. The problem has been discussed in the recent 
literature []6|,[T0HT^1 • The prescriptions offered reduce to a set of rules which are applied to 
the integrands, before the integration is carried out. We illustrate it in the following section. 



V. THE TWO-POINT FUNCTION 

In order to illustrate the method that allows one to subtract the uninteresting polynomial 
part, we consider here in detail the two-point function of the heavy field at order e 2 , in the 
absence of the external field. 



A. The two-point function at p=0 

We first consider again the case where the three momentum is set to zero, and start with 
the expression ( 4.1 5| ) . Combining the denominators by use of the Feynman-Schwinger type 



formula 

1 



ab Jo (a + by) 2 



one may perform all integrations, with the result 



e 2 



S ^ = "l6VM {$llna + $2} + OP " 4) ' 
$i = l--^, $ 2 = 2[87r 2 L A / + ln2-l + ^ln(l + iy)] , W = . (5.2) 

In order to remove the pole in $2, a mass counterterm is needed. In order to keep the 
position of the pole at p° = M, the finite part of the counterterm has to be appropriately 
tuned. In addition, it is explicitly seen that the self-energy develops a logarithmic branch 
point at p° = M. This singularity is suppressed, in the sense that it is proportional to a. 

We now compare this result with what one obtains through the multipole expansion 
PJTH|-p^|. In the present case, this method amounts to expand the integrand in ( 4.15| ) in 



the three-momentum of the massless particle, according to the counting 

p° - M = 0(v 2 ), \p\ = 0(v) , 1° = 0(v 2 ), \i[ = 0{v 2 ) , (5.3) 

where v is a small parameter. Proceeding in this manner, an infinite series of terms is 
generated out of the quantity 

00 

n=0 

e 2 1 rdPl ( l 2 V 1 



Sn ^ ^ AM 2 i(2tt) d J -Z 2 UmJ (M-p° + l°) n + 1 ' ^ 
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The hat on indicates that the result of this manipulation is not identical to the previous 
expression T l f fR . It is seen that £„ is of order v 2D + 2n ~ e near threshold. Let us consider the 
first term in the series, 

e 2 r(I>-2)r(3 n, D _ 3 



4M 2 (47r)^r( 



(M-pT ■ (5.5) 



In contrast to the case considered before, this term is now suppressed near threshold, being 
of order v 2 at D — > 4. The mass is unaffected by the loop correction, and the counterterm 
needed to cancel the ultraviolet divergence is also of order v 2 , provided that we count the 
differential operator idt — M + ^ as a quantity of order v 2 . Furthermore, the non-analytic 
part in this first term, 
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alia a + ■ ■ ■ , (5.6) 



exactly agrees with the non-analytic part in the original expression ( |5.2j ). This feature is 
true to all orders in the v 2 expansion ( |5.4J) . Indeed, the integrations can be done in each 
term, by integrating first over 1° in the lower half plane. In the limit D — > 4, the summation 
can be explicitly performed, with the result 

fii = $i, n 2 = (l-w) {87r 2 L M + ln2-l-ln(l + W0} . (5.7) 

It is seen that indeed, the logarithmic part in T,^ R is exactly reproduced, whereas the second 
contribution <3> 2 , analytic at threshold, is modified, as a result of which E n is of the order 
expected from counting powers of v in the integrand. [The summation of the series in 
Eq. (|5.4| ) has been carried out for illustrative purpose only - in order to demonstrate that 
the infrared-singular part of the multipole-expanded self-energy diagram indeed reproduces 
the exact result given by Eq. (|5.2| ) in all orders in the expansion in a. In practice such a 
summation is never used - instead, only a finite number of terms is retained that contribute 
to the matching condition in a given accuracy] 

The self-energy ( p.7| ) has, aside from the branch point at p° = M, in addition a branch 
point at p° = M/2 on the first sheet in the complex p° plane. This additional branch point, 
however, lies outside the region p° ~ M where both the non-relativistic approach and the 
multipole expansion are supposed to be meaningful. The close analogy with the "infrared 



regularization" of Feynman integrals in Baryon ChPT |15[ is evident. Further, it is seen 
from Eq. ( |5.7|) that the ultraviolet-divergent term is no longer polynomial in momenta after 
the summation - the renormalization procedure can be carried out only order by order in 
the v 2 expansion. 

We conclude that in this example, the two expressions and differ by a power 
series in a. Next, we investigate the contributions from the graphs Fig. [5|B + |5]C. Proceeding 
analogously, it is seen that these start to contribute now at order <j D ~ l . We have not carried 
out the summation of the multipole expanded terms in this case. However, we have checked 
that the leading term in the logarithmic singularity again agrees with the one in the original 
expression in Eq. ( [4.161) . 
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We see that the multipole expansion eliminates the unwanted terms that do not respect 
power counting. In addition, it is clear from the procedure that contributions generated by 
the higher order terms in the effective Lagrangian will be suppressed in the momentum ex- 
pansion in an analogous manner. At each order in the expansion, the ultraviolet divergences 
may be removed by adding a local term to the effective Lagrangian. These terms again 
are suppressed by inverse powers of M, as a result of which one has a perfectly consistent 
procedure. The self-energy T,nr becomes in the multipole expansion 

oo 

Z NR = £ a n {^ln a + d n } + 0(D - 4) . (5.8) 

n=l 

The terms of order a n are generated by a finite number of diagrams. Further, the branch 
point in the propagator stays at M at each order in this expansion, there is no need to 
adjust the mass term in the Lagrangian for this to achieve. 



B. The two-point function at p ^ 0. 

So far, we have considered the propagator at zero three momentum. Let us now consider 
the general case, 

G NR {p) = i J d D xe ipx (0\Th{x)rf(0)\0) . (5.9) 

As before, we consider the corresponding self-energy at order e 2 . The complication that 
arises for non-vanishing three momentum is the following. The contributions from the mass 
insertions h^(A/M 2 ) n h generate singular contributions near p° = M + p 2 /(2M) for small 
three momenta (see Fig. [(]), of the form 

D 2 

const x e 2 p 2 7(M + ^- - p°) m . (5.10) 

The origin of this singularity becomes very clear in the multipole expansion. Indeed, sum- 
ming up the most singular parts of these divergences, one finds that the self-energy behaves 
as 

e 2 

W^P 2 ) = -[Q^M ZlnZ+ "' ' z = (u p -P°)/M , (5.11) 

near the threshold z = 0. Whereas the lowest order diagram Fig. ||A generates a branch 
point at p° = M + p 2 /(2M), the mass insertions move this to p° = u p , as it is required from 
the fact that the theory should describe the non-relativistic limit of the relativistic theory 
pl| . This is analogous to the free propagator, that has, in the absence of any interaction, a 



pole at p° = M + p 2 /2M - only after an infinite re-summation of mass insertions does this 
change to the required form l/(u p —p°). 

The calculation of the non-relativistic propagator may now be performed as follows. All 
calculations are done using the multipole expansion, counting M + p 2 /(2M) — p° and p 2 as 
quantities of order v 2 . At each order in this expansion, the counterterms required may be 
chosen such that Eat_r remains finite at D — > 4, 
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W^,p 2 ) = t ^E^K^- 4 + <} 



16vr 2 M 



e 



n=l 

2 



167T 2 M 



^ n {c n lnz + d n }+0(L>-4) , (5.12) 



n=l 



where the coefficients c' n , d' n and c n , d n are functions of the three- momentum, c' n = c' n (p 2 ), 
etc. This procedure generates a well defined series at each order in the expansion, and each 
term in the series (|5.12| ) is again generated by a finite number of diagrams. The behavior of 
the two-point function near the pole at p° = uj p is given by 

Gnr(p) = Q I ; - ^\ + 0(z D - 5 ) , (5.13) 

u p - p° - E NR (z, p 2 ) ^ P - P° 



with 



^ = l + T ^rf' 1 + 0(e 4 ). (5.14) 



Equation ( 5.13 ) may be used to define the mass in the non-relativistic theory - it is the 
parameter M with the property that the two-point function has a pole at p° = y/M' 2 + p 2 . 



For a detailed discussion of this point, see Ref. [21 



For later use, we give the explicit expression of the multipole-expanded self-energy in the 
non-relativistic theory at order e 2 t> 2 . At this order, only the graph shown in Fig. [5]A con- 
tributes. One should add the counterterm that is needed to cancel the ultraviolet divergence 
- the corresponding Lagrangian is given by Eq. (|6.1|) . The self-energy is given by 

~ e 2 (Mz) D - 3 T(D-2)T(3-D) e 2 f lZ 

NR ~4M* ( 47 r)V T(^±) AM 

= 16^M { ln 2Z + ln 7 - 1 ~ fM l + ° {D ~ 4) ' (5 - 15) 
where at 0(v 2 ), Mz = M + p 2 /(2M) - p°, and 

fi = ^- 2 CL^) + fM) + 0(e 2 ). (5.16) 

From Eq. ( |5.15| ), one can read off the expression for the wave function renormalization 
factor, 

^ = l-0 + O(eV,e 4 ). (5.17) 



C. Comparison with the relativistic theory 

It is instructive to compare at this stage these results with the two-point function of the 
relativistic theory considered in section ||. In the absence of the external field, one has 
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i J Ae- (0\TH(x)W(0)\0) = w _ J_ s(p2) ■ (5.18) 

The relevant contributions at order e 2 are displayed in Fig. |l]A, |lp. Adding the counterterms, 
we find 

S(p 2 ) = j^j- 2 r(2 - £) jf 1 rfx [xM 2 - x(l - x)p 2 ] D / 2 - 2 - 5M 2 . (5.19) 

Performing the limit D —>■ 4 before going to the mass-shell, we obtain 

^,^^ln^ + (.-4,. (5,0, 

Now, changing the sequence of the limiting procedures and differentiating E at / = M 2 
and .D > 4, one obtains the wave function renormalization constant Z, 

~ ' ' L M + -^-- + 0^). (5.21) 



2M 2 M 16vr 2 M 2 

By comparing the leading logarithmic singularity in S with the one in given by 

Eq. (|5.11|) , it is seen that they are identical, up to a factor l/2co> p . Based on the results in 
this section, we expect that higher order non-analytic terms in the relativistic self-energy 
will be reproduced by the multipole expansion as well. The polynomial part can then be 
adapted by a suitable choice of the counterterms, as a result of which one obtains 

± NR (z,p 2 ) = ^- S(p 2 ) , p 2 ~ M 2 , (5.22) 

order by order in the v expansion. We refer the interested reader to Ref. for a discussion 
of this relation beyond one-loop accuracy. 



VI. POWER COUNTING AND MATCHING 

We now introduce the external field <fi again, and consider radiative corrections of order 
e 2 (one loop) to the scattering matrix elements discussed in section III. We restrict the 
discussion to Green functions that are relevant for bound state calculations considered later 
on - i.e., we consider the one heavy particle sector in Fock space, with no incoming or 
outgoing light particles. 



A. The Lagrangian, the reduction formula and power counting 

To generate the tree graphs, one replaces the external field A0 in the Lagrangian given 
in Eq. (|C2j), Xcf) — > X(f> + e£, and adds the kinetic term for the light field. Because we do not 
consider incoming/outgoing light particles, one needs to retain in the interaction Lagrangian 
only terms that are linear in the field I. In the presence of radiative corrections, as we have 
already seen, one has to equip this Lagrangian with counterterms that render the theory 
ultraviolet finite, 
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/ A \ 1 00 1 

C m = ft* ( id, - M + — ) h + + £ ^ £„ + £ , , 

C 3 = ft*(A 2 + 93 A(^A + A0) + 94 AV + <j 5 e(«A + AC) + 2g e e\<j>t)h , 

£ "=SP ,I, ( '- M + lf> + "'- < 61) 

The coefficients gi appearing in the Lagrangian are power series in the variable e 2 /M 2 , 

Qi = 1 + 0(e 2 ), and f\ was determined in (|5.16|) at leading order in e 2 . Scattering matrix 

elements are calculated by use of the reduction formula Eq. (|3.8| ) except that now - in analogy 

to the relativistic framework - the external legs of the heavy particles are equipped with wave 

1/2 

function renormalization factors Z^ R , and the radiative corrections to the external legs are 
omitted. The matrix elements are then expanded in powers of the coupling constants A and 
e, and of the external momenta. We define the power counting scheme - that specifies how 
the momentum expansion is done - by associating the following powers of v to the derivatives 

idt-M + ^ = 0(v 2 ) , A = 0(v 2 ). (6.2) 

[Note that one may arrange the derivatives in such a manner that the Laplacian never acts 
on the light field.] In addition, it is understood that the loops are evaluated by use of the 
multipole expansion considered above: the loop momenta of the light particle are counted 
as quantities of order v 2 , and the propagators of the heavy field are expanded accordingly. 

One could now set up a formal power counting for generic graphs containing generic 
vertices from ( |6.1|) . However, we find it more instructive to consider instead two examples, 
the self-energy discussed in the previous section, and radiative corrections to scattering in 
the external field at lowest order. 

We start the discussion with the self-energy - a few graphs are displayed in Fig. 
According to (|6.2|), we count the heavy (light) propagator as a quantity of order l/t> 2 (1/f 4 ), 
as a result of which the graph Fig. §A is counted as a term of order v 2 . As we have discussed 
in the previous section, the multipole expansion transforms the relevant integral into an 
infinite string of terms, starting at 0(v 2 ). In addition, this graph requires, in the multipole 
expansion, a string of counterterms that start with the one proportional to f\ in ( |6.1| ). This 
is indeed a term of order v 2 according to the above counting. The vertices in Fig. [5|B, ||]C stem 
from £3 and generate therefore additional factors of at least the order v 2 . The corresponding 
string of counterterms starts at order v 4 , etc. We now turn to scattering in the external field 
at order A. Some of the relevant graphs are displayed in Fig. |8|. The tree graph Fig. has 
been evaluated in ( |3.12| ). The radiative corrections in Figs. §B and |8p amount to multiply 



the tree graph with the residue Z NR of the propagator - this contribution starts at order 
e 2 \v°. The vertex correction Fig. |8p starts at the same order (we discuss below the manner 
in which one tames the ultraviolet and infrared divergences hidden here) . Finally, the graphs 
Fig. ^E and Fig. |8|F contain one vertex from £3, proportional to g§. The leading term of 
this contribution is suppressed by v 2 , because it contains one heavy propagator only. 
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These two examples show two things: First, one associates to any tree or one- loop graph 
a factor v m according to the rules spelled out above. The multipole expansion turns the 
loop contribution into an infinite string of terms, such that the amplitude starts at 0(v m ). 
Second, contributions from higher order operators C n are suppressed in the following sense. 
Consider loop corrections to a matrix element at some fixed order A n , and compare the 
corrections due to a loop generated by C\ x C p with the one generated by C\ x C q , with 
p > q. Then, the leading contribution from the former is suppressed with respect to the one 
of the latter, and the relevant string of counterterms needed the renormalize the former starts 
at higher order than the one for the latter: although individual graphs do not contribute at 
a single order, the starting point of the expansion is hierarchically ordered. 



B. Matching: The vertex function at order e 2 

In order to determine the coefficients gt in flS.ip, one evaluates the amplitudes at a fixed 
order in A, e and requires that the results be the same in the relativistic and in the non- 
relativistic theory at a given order in v. As the effect of the operators of higher dimensions 
are suppressed in the multipole expansion, only a finite number of coefficients contributes 
at a fixed order v. In the following section, we will show that we need only the coefficient 
gi at order e 2 for our purpose. It suffices to evaluate a matrix element that fixes this term, 
and we choose scattering of the heavy particle in the external field <fi at order e 2 A. 

The relevant diagrams in the relativistic theory are displayed in Fig. [7| (with counterterms 
omitted), while those of the non-relativistic scattering matrix T^r are shown in Fig. || As 
we are considering a massless theory, we use dimensional regularization to also regularize 
the infrared divergences. 

The matching condition of relativistic and non-relativistic scattering amplitudes reads 

7V*(p,k)=T(p,k). (6.3) 

At order A, these scattering amplitudes can be written as 

T(p,k) = A0(p-k)T(p,k), 
T NR (p, k) = A0(p - k)f NR (p, k) . (6.4) 

At D ^ 4, both T and Tnr can be expanded in the momenta. Since we only need to 
determine the coefficient g±, it suffices to work at vanishing momenta, and to require that 

f NR {0, 0) = f (0, 0) . (6.5) 

We start the calculation with the relativistic amplitude. The contribution from diagram 
Fig. [7|E has been evaluated in section |J. The diagram Fig. |7p is ultraviolet finite and easy 
to calculate at zero three-momenta in D dimensions. Adding the contribution from the wave 
function renormalization constant Z gives 
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As expected, the infrared divergences that are present in the residue Z and in the vertex 
function Fig. [7|D, cancel at threshold. Indeed, infrared finite cross sections are obtained by 
adding diagrams from Bremsstrahlung. At threshold, the phase space for these processes 
vanishes - therefore, the amplitude must be infrared finite itself. 

Next we calculate the amplitude in the non-relativistic theory. The vacuum polarization 
diagram is absent in this case, while the contributions from Fig. |8]E, |F are suppressed 
in the v expansion. The contribution from diagram Fig. [5p is ultraviolet divergent in the 
multipole expansion. In order to separate this divergence from the one in the infrared region, 
we proceed as follows. We first consider the amplitude off-shell, 

T NR (p°) = 4M 2j( 27r )D / _p (m - p° + i°y ' ^ 
Evaluating this expression in D dimensions gives 

e 2 rg-2)T(4-j) 0D _ 4 

Now, we eliminate its divergence at D = 4 by using 

9i = l + (£(/,) + fM) + 0(e 4 ) , (6.9) 

where g[{fi) denotes the finite part in the coupling constant g\ at order e 2 . Finally, keep- 
ing D > 4, multiplying with Z^r the lowest-order contribution displayed in Fig. ||A and 
summing up all contributions, we obtain 

f NR (0, 0) = 1 + (g r M - fM) + 0(e 4 ) . (6.10) 

[Remark: under the following sequence of limiting procedures: going to mass- shell p° -> M 
at D > 4 before performing the limit D — > 4, the unrenormalized vertex function Tjv_r 
vanishes. It is only the polynomial part proportional to g%, that survives in this limit.] 

The matching condition ( |6.3| ) now determines the difference of the finite parts of the 
counterterms, 

9 r M-f r M = l + 0(e 2 ). (6.11) 



C. The Hamiltonian 

In order to calculate the energy-levels, we need to construct the Hamiltonian pertaining 
to the Lagrangian (|6.1|). As is discussed in Appendix 0, we first re-scale the heavy field, 
such that the kinetic term has coefficient one. As we already mentioned and as is proven in 
the following section, only the terms proportional to gip are needed in the calculation of the 
energy-level at order e 2 A 2 . Further, terms of order e 2 in gi do not contribute. Therefore, 
re-scaling the heavy field amounts to the renormalization 
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e 2 5e 2 
9l ^ 9l -UP fl = 1 + 96^ ' (6 ' 12) 



result of which the Lagrangian becomes 
C NR = h^id t -M+^h+ l - (d,£) 2 + 



A / 5e 2 

+ 777 1 + 



2M 2M V 96vr 2 M 2 



(6.13) 



The ellipsis denotes terms that are not needed in the following. 

Remark: This Lagrangian was obtained without specifying the couplings f\,gi - only the 
combination fl6.12|) counts, which is fixed by the matching condition ( |6.11|) . The Lagrangian 



( |6.13| ) amounts to using fi = 0, which leads to Znr = 1 at this order in the low-energy 
expansion. As can be seen from the explicit expression ( |5.15| ), the derivative of the non- 



renormalized self-energy indeed vanishes at D > 4, which is a particular manifestation 
of so-called "no-scale" arguments. In this setting, the self-energy is not finite at D — 4. 
However, the pole position of the two-point function is not affected, and S'-matrix elements 
remain the same.Q 

Finally, we display the Hamiltonian. Denoting by ILj the canonically conjugated mo- 
mentum of the light field, one has 



with 



n = n + n c + n I , (6.14) 



H ° = h i M -ii) h+1 2 U * +1 2^) 



Uc = ^ h \ht h Hl = h l£ h h ^ h . ( 6 .i5) 

° 2M ' 2M 192tt 2 M 3 y v ; 

We use in the last section these expressions to calculate the energy-level shift of the ground 
state. 



VII. THE GROUND STATE 



The radiative shift of the energy-levels is calculated by use of Eq. ( |3 ■ 3 1| ) 



AE - / |H I \ \ - (w|Hi|m)(m|Hi|ra) y> (w|Hi|m, 7) (m, 7|Hi|ra) 

n yl\ l|"/ ~T / , rp rp ' / 1 rp rp rp 

m -L n - C/ n - C/ m m,7 • C/ n • C/ m - C/ 7 



rn=£n 



(E n — E„ 



4 We are under the impression that many of the contemporary bound-state calculations use im- 
plicitly an argument of this kind. On the other hand, we are not aware that the relation of this 
procedure to standard renormalization theory has been made explicit in the literature. 
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where 7 denotes the light particle which is present in the sum over intermediate states 
in the Fock space. The Hamilton operator is constructed from the Lagrangian displayed 
in Eqs. ( |6.1| ). The interaction Hamiltonian is no longer diagonal and induces transitions 
between the Fock space vectors with an arbitrary number of light particles/antiparticles. 
The unperturbed wave functions in this expansion are the eigenfunctions of Ho + He, 

(H + H c )|n> =E n \n), (H + H c )|n, 7) = (E n + EJ\n, 7) • (7.2) 

The Hamilton operator contains an infinite string of operators. Only a few of them, with the 
low mass dimension, contribute to the energy shift in a given accuracy. For the lowest-order 
term in ( |7.1| ), or for the sums with no light particle in the intermediate state, one can apply 
power-counting arguments: Since the Coulomb wave functions do not depend on the heavy 
mass M, and the energy denominators are homogeneous functions of M, one can count the 
powers of M in the energy shift once one knows the explicit powers of M in the couplings of 
various operators in Hj. Given more powers of M in the denominator, one needs additional 
powers of A in the numerator, in order to compensate the extra mass dimensions - that is, 
only a few first terms in the Hamiltonian can contribute to the energy at 0(e 2 A 2 ). The 
same argument with a slight modification applies also to the sum with one light particle 
in the intermediate state. Summing up the Coulomb wave functions, one ends up with the 
Schwinger's Green function, with possible insertions of higher-dimensional operators, folded 
by the light particle loop. Further, one expands the Schwinger's function in powers of the 
external field - then, each integral has the form already considered above, in the context of 
the scattering problem. Using the multipole expansion restores the power counting - the 
leading term is determined by the naive power-counting rules, and the sub-leading terms 
are explicitly given. Consequently, for this term in the perturbation expansion one can also 
apply the same power-counting arguments that were used in the terms with no light field. 
Therefore, one has to examine only few contributions from the lowest-order operators in Hi, 
in order to obtain the energy-level shift in a given accuracy. In particular, the coefficient 
gi is needed in the accuracy 0(e 2 ), whereas it suffices to determine #2 at tree level, and 
all higher-order terms can be ignored. To conclude, in order to evaluate the energy-level 
shift at 0(e 2 A 2 ), it suffices to work with the Hamiltonian given by Eqs. (|6.14[ )- (|5.15 ). For 



a more detailed discussion of the power counting in the bound-state calculations within the 
framework of non-relativistic effective theories, we refer the interested reader to Refs. ||[7|]. 

Let us now pass to the calculation of the ground-state energy-level shift. The matrix 
element (l|Hi|l) receives contribution only from the second term in the interaction Hamilton 
operator, 



The second and fourth terms in the perturbative expansion ( |7. 1| ) do not contribute to 
the bound-state energy at 0(e 2 A 2 ). In the third term, the summation over the Coulomb 
wave functions can be carried out explicitly, resulting in Schwinger's Green function 



Then, the contribution from this term to the energy-level shift can be written as follows 

Jd» d (p)dv d (lL)tf(j>ytf(k) J ^^(27r) d G s (E 1 -M-|q|;p-q,k-q). (7.4) 



e 2 



AM 2 
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An important remark is in order. In Eq. ( [7.4|) , one has to use the Schwinger's function 
in d dimensions, defined through the spectral sum over (^-dimensional wave functions (see 
Appendix [TJ). Further, this function can be written as a sum of zero-Coulomb, one-Coulomb 
and multi-Coulomb parts. In the first two, the generalization to d dimensions does not cause 
any problems, while the multi-Coulomb part is ultraviolet convergent and can be considered 
at d — 3. 

The zero-Coulomb contribution to the bound state energy is given by 

AE II0 = ( dvJp) |^f(p)| 2 / dvAq) —- -, ^ . (7.5) 

' AM 2 J j d ™ 2\q\ Ef - M - \q\ - 

In order to be consistent with matching in the scattering sector, one has to carry out the 
multipole expansion in this integral as well. The result at 0(e 2 A 2 ) then is 

e 2 r f , , m ,d, ms P 2 + 2M(M - Ef) e 2 A 2 / A 2 3\ . , 

The one-Coulomb contribution to the bound-state energy is calculated analogously, 

A-E/i.i = J di/ d (p) dv d (k) <(p)*<(k) j-l^ 2 x (7.7) 

x f dvJq) —, — r ? \a r, vr- ■ 

J W 2|q| (Ef - M - |q| - ^)(Ef - M - |q| - ^) 

After carrying out the multipole expansion, we obtain 

(7.8) 

Note that the divergent terms as well as terms non-analytic in A, cancel in the sum AEjj^ + 

In the multi-Coulomb part of the energy shift, as it was mentioned above, one can 
directly perform the limit d — *■ 3, since this contribution is ultraviolet and infrared finite. 
After re-scaling the integration variables, this piece is written as follows 

87r 2 e 2 A 2 f 1 1 , , 

AE II;m = j^- J di/(u)dy(v) (1 + u2)2 (1 + v2)2 x (7-9) 

x / du(w)- ' 99 



where 



wu/i + iwi -o F(p;|w|;u,v) 1 + |w| + (u - ^w) 2 1 + |w| + (v - ^w) 



F(p; |w|; u, v) = (u - v) p + ^jrnr (l + |w| + (u - t^w) )(1 + |w| + (v - — w 



(7.10) 
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and 7 = A/2. Taking the limit 7^0, one identifies this expression with the analogous con- 
tribution in the relativistic case, Eqs. ( [2. 41 ) and (p. 42 ). Finally, collecting all contributions 
to the radiative energy-level shift, we obtain 

e 2 A 2 



AE l = (l|Hi|l> + AE IIfi + AE n>1 + AE n , r 



.( 



16 



81n2 + J + ■■ 



64vr 2 M 3 V 3 

which is exactly the result ( |2.43| ) obtained by using the relativistic framework. 



(7.11) 



VIII. SUMMARY AND CONCLUSIONS 

i) We have considered in this article the energy-levels of a massive scalar particle, which 
moves in an external Coulomb field and interacts in addition with a massless scalar 
particle through a Yukawa interaction. This has allowed us to study in detail the non- 
relativistic effective Lagrangian treatment |||] of the bound-state problem in Quantum 
Field Theory. 

ii) We first considered the relativistic calculation of the energy levels. The calculation 
serves as a reference point for the non-relativistic calculation performed later on - with 
the exception of the non-renormalizable interaction considered at the end of section II. 
Indeed, that interaction mimics the situation encountered in ChPT. We have shown 
that - at lowest nontrivial order in the relevant coupling constant g - bound-state 
observables are ultraviolet finite, once the Green functions in the scattering sector 
are made finite by an appropriate choice of the counterterms. Although we have not 
proven it, we have no doubt that this will remain true to all orders on g. 

iii) A systematic framework for the perturbative calculation of bound-state observables 



was then developed in the non-relativistic sector, based on Feshbach's |17j technique. 
At the end, one recovers the conventional Rayleigh-Schrodinger perturbative expansion 
of the energy-levels. The method allows one in addition to calculate the position of 
the pole on the second Riemann sheet, in case that the bound state turns out to be 
unstable. 

iv) Although the external field problem may sound trivial at first, actual calculations in the 
non-relativistic theory reveal that this is not the case. Indeed, applying the Rayleigh- 
Schrodinger perturbative expansion, one has to verify that the ultraviolet divergences 
in the matrix elements of high-dimensional operators cancel, provided that all con- 
tributions at a given order in coupling constant are summed up. Stated differently, 
the non-relativistic expansion does not generate additional ultraviolet divergences, de- 
spite the fact that the non-relativistic Lagrangian contains terms with arbitrarily high 
powers of spatial derivatives. 

v) We have explicitly checked in an example that the off-shell ambiguity in the effec- 
tive theory - that can be traced to the presence of terms in the Lagrangian which 
vanish upon use of the equation of motion or can be removed by field redefinitions 
- does not affect the energy-level. [Similar results were recently obtained in the con- 
text of NRQCD [p2] and in the non-relativistic effective field theory approach to the 



deuteron 
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vi) Next, we considered radiative corrections in the non-relativistic sector in some detail. 
The construction of the relevant effective Lagrangian is based on the following idea. 
First, one constructs a Lagrangian that agrees at tree level with the relativistic ampli- 
tude to all orders in the low-energy expansion. This guarantees that the correct cut 
structure is obtained at one-loop level, in the sense that the absorptive parts of the 
relativistic theory are reproduced. Loops then simply amount to close two of the ex- 
ternal lines in the tree graphs - therefore, we expect that the non-relativistic ^-matrix 
elements do have the correct non-analytic structure at one-loop order. 

vii) We have investigated the two-point function of the heavy particle at second order 
in the Yukawa coupling. It was demonstrated that - at zero three-momentum - the 
multipole expansion of the Feynman integral reproduces the non-analytic parts of 
the original expression, order by order in the expansion. The polynomial part of the 
integral - which leads to the breakdown of counting rules - is modified by the multipole 
expansion. This merely amounts to a change of the renormalization prescription and 
does not affect physical observables. 

viii) In addition, we have investigated the structure of the two-point function in the thresh- 
old region, and have shown how the correct pole structure emerges from the multipole 
expansion. Eq. Q5.22| ) displays the expected relation between the non-relativistic and 
the relativistic self-energies at one-loop order, but at any order in the f-expansion. 
That relation presumes that one has properly chosen the counterterms. 

ix) We have considered the renormalization procedure for the non-relativistic effective 
theory, including the subtleties that arise due to the presence of zero-mass particles. In 
particular, the renormalization of two and three-point functions in the non-relativistic 
theory was treated in detail. Further, we have demonstrated that - at leading order 
in the low-energy expansion - a procedure which utilizes no-scale arguments in the 
dimensional regularization is equivalent to the standard renormalization scheme. 

x) Constructing finally the non-relativistic Hamiltonian in the presence of the light field, 
we have verified that it reproduces - at first nontrivial order in the Yukawa coupling - 
the relativistic expression for the energy-level of the ground state. 



In conclusion, we find that dimensional regularization, together with the multipole ex- 
pansion []6|,p~Q| p~2|| and Feshbach's representation [17| of the Green functions, allows a sound 
evaluation of energy-levels in the non-relativistic models considered. Furthermore, with the 
same technique, one can determine the position of the poles on the second Riemann sheet 
for decaying states Pfl. No use was made of a representation of the Lagrangian in terms 



of expanded fields |]10|,|T2|] in order to perform these calculations. 



Acknowledgments 

We are grateful to M. Beneke, H. Griesshammer, M.A. Ivanov, T. Kinoshita, P. La- 
belle, G.P. Lepage, H. Leutwyler, V.E. Lyubovitskij, K. Melnikov, P. Minkowski, M. Nio, 
H. Sazdjian and J. Soto for useful discussions. This work was supported in part by the Swiss 



35 



National Science Foundation, and by TMR, BBW-Contract No. 97.0131 and EC-Contract 
No. ERBFMRX-CT980169 (EURODA$NE). 



APPENDIX A: ONE-COULOMB CONTRIBUTION TO THE ENERGY-LEVEL 

SHIFT 

In order to calculate the one- Coulomb contribution to the energy shift, we first combine 
the denominators in Eq. ( |2.39| ) by using Feynman parameters. Then, one obtains 

1 f l In g(x) — In f(x) 



r(E 1 ;p,k) = -— / dx 

l07T z JO 



167r 2 Jo g(x) — f(x) 

g(x) = M 2 + x(l-x)(p-k) 2 , f(x) =x(7 2 + p 2 ) + (l-x)( 7 2 + k 2 ). (Al) 



Next, we substitute Eq. (|A1|) into Eq. (|2.38|) , and re-scale the integration variables ac- 
cording to p = 711, k = 7V. We obtain 

F n>1 = 8z 2 (ln z 2 - 2 In 2 - 1) h + 8z 2 I 2 + o(z 2 ) , (A2) 

where 

h = J dv(u)dv(v) ^^- 2 = JL- , (A3) 

and 

1 (l + u 2 )ln(l + u 2 ) - (1 + v 2 )ln(l +v 2 ) 



J a = / dv(u)dv(v) — l — - 



+ u 2 ) 2 (1 + v 2 ) 2 |u - v| 2 u 2 -v 2 

(A4) 

In order to calculate I2, it is convenient to use the integral representation 
(1 + u 2 ) ln(l + u 2 ) - (1 + v 2 ) ln(l + v 2 



Jo \x + l 



u 2 -v 2 Jo \x + l (x + l + u 2 )(x + l + v 2 )/ ' 

(A5) 



and to introduce the Fourier transform of |u — v| 2 , 



1 1 r p «(u-v)r 

J_ \ d ^ . (A6) 



|ll — v| 2 47T 

After this, all integrations in I 2 can be performed analytically, resulting in 

h = -\n2] . (A7) 

128tt 2 V2 J K J 

Substituting Eqs. ( [A3] ) and ( |A7| ) into Eq. (|A2|) , we arrive at the final result displayed in 
Eq. (F40|) . 
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APPENDIX B: MULTI-COULOMB CONTRIBUTIONS TO THE 

ENERGY-LEVEL SHIFT 



In this Appendix, we present the calculation of the multi-Coulomb contribution to the 
energy-level shift. In this calculation, a straightforward use of Schwinger's representation en- 
counters a difficulty, since in the integration region — oo < q° < oo the variable v in Eq. (|2.14|) 



can become larger than 1. In order to overcome this difficulty, we use the spectral repre- 
sentation for the Schwinger's function in terms of the solutions \n) d of the non-relativistic 
Schrodinger equation with the Coulomb potential (|3.34j ). By using this representation, the 
integration over the variable q° can be performed with the use of Cauchy's theorem. The 
whole contribution from the diagrams of the type II to the self-energy part in D dimensions 
can be written in the following form 

r n (F d -r, lc\ = -,p 2 V f d ° q 1 (P-ql w )<* d{n\k-q) 
1 1,P ' >~ (2ir) D q 2 + tO (Ef-q°) 2 -(E d ) 2 + tO 

+ EiTWTU> , (B1) 

where the sum runs both over the discrete and the continuous spectra of the unperturbed 
Coulomb Hamilton operator. 

The representation (|B1|) now includes all contributions from the diagrams of the type 



(P - q|«)d d(n 


|k-q) 


( 1 


AE d 

n 


q| 




^E*-Ef + 


q| 



II, apart from the mass counterterm — 5M 2 . To single out from this expression the low- 
energy contribution which contributes at the leading order in A, we introduce an additional 
integration over the auxiliary variable x, using the following identities 



1 


29 (Ef - 


q| 




f°° dx (Ef - 


q|) 


E n( E n ~ Ef + |q|) (E*) 2 -(Ef 




q|) 2 it J 


o (x 2 + (Ef - 


q|) 2 )(* 2 + (^) 2 ) 



dx (Ef+ |q|) 



E*(E* + Ef+\q\) nJo (x 2 + (Ef + |q|) 2 )(x 2 + (E%) 



2\ 



(B2) 



If now one substitutes Eq. (|B2j ) into Eq. ([BID, one can single out the piece corresponding 
the multi-Coulomb exchanges. [In this piece, one can safely put D = 4, Ef = E u etc, since 
there are no ultraviolet divergences.] As far as this contribution is concerned, it is seen that 
only the part of this expression which does not involve the additional integration over x (the 
one coming from the first term in the r.h.s. of the first line of Eq. (p2|)), contributes at order 
0(X 2 ). Other terms can be proven to contribute at order 0(A 4 ). The reason for the use of 
the above trick is the following. One can observe from Eq. ( |B2| ) that now the variables which 
serve as an "energy variable" in the Schwinger's representation: e = ((E± — |q|) 2 — M 2 ) / (2M) 
for the term without x integration, and e(x) = —(x 2 + M 2 )/(2M) for the terms with the 
x integration, vary in the regions —M/2 < e < — A 2 /(8M) and — oo < e(x) < —M/2, 
respectively. The corresponding values of the parameter v in the Schwinger's representation 
never exceed 1 and, consequently, the integral over dp in Eq. (|2.14j ) always converges. 



According to the above discussion, with the use of the Schwinger's representation one can 
now write down the part of the multi-Coulomb contribution to the heavy particle self-energy 
graph, which gives the leading-order contribution to the energy shift 
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£ (4p-q,k-q)-ej (M 2 — (Ei — Iq]) 2 ) 1 / 2 X 

X M 2 + (p — q) 2 - (E l - |q|) 2 M 2 + (k — q) 2 - (E 1 - |q|) 2 + " ' ' (B3) 

where the quantity I is given by Eq. ( |2.14j ), and the ellipsis stands for terms which generate 
higher-order contributions in A. 

Substituting Eq. ( |B3[ ) into the expression of the energy-level shift and re-scaling the 
integration variables p = 711, k = 7V, q = 7 2 w/(2M), one arrives at the Eqs. ( |2.41| ) and 

(ED- 
APPENDIX C: USE OF THE EQUATION OF MOTION 
1. The effective Lagrangian ( |3.6[ ) 



The original form of the Lagrangian is obtained by expanding the propagator Dz 1 and 
the differential operator d - that appear in the operator T> in ( |3.2j ) - in inverse powers of the 
heavy mass, 

^--{^M+v^^^j^^^K 2 ^)^.], 

rf = ^=( 1 + 7^ + ---l ; d t = id t -M. (ci) 



V2M I 4M 2 
The result is 

C NR = tfvh = h^D + h + £ pr^ (n) , 

= Xh^h , £ (2) = , 
£( 3 ) = \h\<j)A + A(f))h + \ 2 h ] <p 2 h , £ (4) = -\ 2 h ] <pd t <ph , 

= \rf (A0A + ^A 2 + ^A 2 ^ h 

+ \ 2 rf (0 2 A + A0 2 + 30A0 + <f)d 2 <p) h + \ 3 rf<f) 3 h. (C2) 

The time derivatives of second and higher order may be eliminated in the following manner. 
First, we write the Lagrangian as 

C NR = tiDh , D = d t + A + A x d t + £ A n d n t . (C3) 

n>2 

The operator A contains Laplacians and external fields 0, and ^1,2,3... contain at least two 
external fields. The two-point function is obtained by evaluating the classical action 

S = J d 4 x(£ NR + rfj+^h) (C4) 



38 



at the solution to the classical equation of motion, 



Dh = -j , (C5) 

result of which 

S = J d 4 xj ] h , h = -D~ x j . (C6) 
We eliminate the time derivatives successively from Eq. ( |C5| ) . First, we consider the quantity 

S = J d 4 xj r h , (C7) 

where 

h = -(^-L^-D)-^ = ~ D ~ lj ~ D ~ lAlJ ■ (C8) 

The second term does not contribute to on-shell matrix elements, because on one of the 
external legs, there is no pole: the source j couples directly to the external fields (j) present 
in Ai. We conclude that the quantity S contains the same on-shell matrix elements as the 
original classical action S. Stated otherwise, one may use, instead of ( |C3|) , the Lagrangian 

C NR = h^D 1 h , (C9) 

with 

Di = — l —D 
1+A 1 

= d t + B +J2B m d? ; B m = -^—A m , m = 0,2,3,.... (CIO) 

m>2 1 

The on-shell matrix elements are obtained by inserting in ( |C6[ ) the solution of 

Dih = -j . (Cll) 

By acting with B^dt on ( |C11|) , solving for B^d\ and proceeding as before, one may eliminate 
B<id\, and similarly for all higher time derivatives. 



2. The effective Lagrangian ( [4.12 ) 

The main observation is the fact that, as far as tree graphs in the one heavy particle sector 
are concerned, we may perform in the classical action (pi) the replacement — > eD^ 1 f, 
where / is the external field associated with the light particle, and □Z)(^ 1 (x) = 5 A (x). 
Therefore, the problem of eliminating higher order time derivatives is reduced to the one 
the previous discussion, except that now, the external field is time-dependent. The time 
derivatives act on the heavy and on the external field. The ones acting on the heavy field 
may be removed as is described above. As for the time derivatives on the external field, one 
uses □D ( ^ 1 (x) = 5 4 (x) and can thus eliminate all of them, except those of first order. 
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3. Re-scaling the heavy field 



Constructing the Hamiltonian, it is useful to rely on a Lagrangian whose kinetic term is 
properly normalized. This is not the case for the one displayed in ( |6.1|) due to the presence 
of the counterterms. Since these are needed only at tree level, one may again use the 
above described manipulations to eliminate them. However, some of the unwanted time 
derivatives differ from the ones in ( p3|) , to the extent that the operators A n need now not 
contain external fields - an example is the counterterm proportional to f\ in ( |6.1| ). As a 
result of this, the manipulations described after Eq. (|C3|) lead to Green functions that differ 
from the original ones also on the mass-shell - to quote an example, the residue of the two- 
point function is be modified. However, one can convince oneself that, in the calculation of 
the S'-matrix elements, this change in the two-point function is exactly compensated by a 
corresponding change of the coupling constants in the effective Lagrangian. As a result of 
this, the manipulations leading from ( |C5| ) to ( |C10| ) generate a Lagrangian that leaves the 
S'-matrix elements untouched. 



APPENDIX D: CONNECTION OF THE RESOLVENT WITH THE TWO-POINT 

FUNCTION 



In this Appendix we derive Eq. (|3.22|) of section III , which gives the connection of the 
resolvent 

with the two-point function 

G NR {E;p,k) = i J dtdVye lEt ~ ipMy (0\Th(t,x)h\0,y)\0). (D2) 

Expressing the creation and annihilation operators of one-particle states in terms of the 
fields, Eq. ( |3.17| ), and using the representation 

E — ti i Jo 

which is valid for Im E > 0, as well as the fact that H generates translations in time 

h(t, x) = e iHt h{0, x)e~ im , (D4) 

we get 

<p|-= — =|k> = - / dt d 3 xd 3 ye lEt - ipMy t(0\e- im h{t,x)h\0,y)\0). (D5) 
E — H Jo J 

We assume that the full Hamilton operator annihilates the vacuum H|0) = 0. Since t > 0, 
we may replace the expectation value of the ordinary product with that of the time-ordered 
product of fields without changing anything. Finally, observing that (0|/i^(0, y) = 0, we can 
extend the t-integration to — oo and find 

(Pl^jjlk) = -G NR (E;p,k). (D6) 

Note that the last step is a consequence of the non-relativistic nature of the theory: the 
number of particles associated with the field h cannot change. 
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FIGURE CAPTIONS 



FIG. |]. The ultraviolet divergent diagrams in the relativistic theory at one loop. The 
heavy (light) particle is denoted by a solid (dashed) line, and the wavy line represents the 
external field. (A) Tadpole diagram, (B) the interaction of the light field with the external 
field through a heavy loop, (C) the light particle self-energy graph, (D) the heavy particle 
self-energy graph. 

FIG. |2]. The self-energy diagrams of the heavy particle moving in an external field (rela- 
tivistic theory, order e 2 ). The counterterm contributions are not not shown. The double 
line corresponds to the propagator G(x,y), and the dashed line denotes a light particle. 

FIG. |3|. Scattering of the heavy particle in the external field at order A 2 . The heavy 
particle is denoted by a solid line, and the wavy line represents the external field. Graph (A) 
contributes to the Green function in the relativistic theory, whereas (B) and (C) contribute 
in the non-relativistic formulation. Crosses denote mass insertions. 

FIG. f|. The scattering process k + q± + ■ ■ ■ + q n — > p + q n+ i + • ■ ■ + q m in the relativistic 
theory at tree level. The heavy (light) particle is denoted by a solid (dashed) line. 

FIG. The self-energy graph of the heavy particle at order e 2 , non-relativistic theory. 
The heavy (light) particle is denoted by a solid (dashed) line. (A) The lowest-order graph 
with non-derivative vertices, (B) and (C) are graphs with one derivative vertex, denoted by 
a filled circle. 

FIG. ^ Mass insertions (indicated by crosses) in the self-energy of the heavy particle, 
non-relativistic theory. The heavy (light) particle is denoted by a solid (dashed) line. 

FIG. [7|. Scattering of the heavy particle in the external field at order e 2 A, relativistic theory. 
The heavy (light) particle is denoted by a solid (dashed) line, and the wavy line represents 
the external field. Diagrams containing counterterms are not shown. 

FIG. ||. Scattering of the heavy particle in the external field at order e 2 A, non-relativistic 
theory. The heavy (light) particle is denoted by a solid (dashed) line, and the wavy line 
represents the external field. Diagrams containing counterterms are not shown. 
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